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Abstract 

It has long been recognized that the dynamics of linear quantum systems is clas- 
sical in the Wigner representation. Yet many conceptually important linear problems 
are typically analyzed using such generally applicable techniques as influence func- 
tionals and Bogoliubov transformations. In this Letter we point out that the classical 
equations of motion provide a simpler and more intuitive formalism for linear quantum 
systems. We examine the important problem of Brownian motion in the independent 
oscillator model, and show that the quantum dynamics is described directly and com- 
pletely by a c-number Langevin equation. We are also able to apply recent insights 
into quantum Brownian motion to show that the classical Fokker-Planck equation is 
always local in time, regardless of the spectral density of the environment. 
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It has long been recognized that, although the Wigner function of a linear quantum 
mechanical system cannot in general be considered as a classical ensemble density, 
yet it evolves according to exactly the same Liouville equation as the corresponding 
classical system: / = L c z/.|1| This is true even for linear systems with time-dependent 
parameters, and for non- linear systems treated in Gaussian approximation 0. Nev- 
ertheless, conceptually important problems, such as quantum Brownian motion and 
particle production in classical background fields, are generally treated using for- 
malisms which do not explicitly take advantage of classicality, such as path integrals 
and Bogoliubov transformations. An important example of the disadvantages of this 
unnecessary sophistication is the way in which the machinery of influence functionals 
has successfully hidden the fact that, in the independent oscillator model |§, quantum 
and classical Brownian motion are dynamically identical. 

The more sophisticated approaches are often technically cumbersome, and their 
interpretation frequently seems rather subtle. It is usually a challenge, for example, 
to associate a simple physical intuition with the closed time paths of an influence 
functional. In contrast, a much more direct approach to any linear problem is to use 
the classical equations of motion to obtain the time evolution of the Wigner function. 
After solving these equations, with their (usually easiest) Cauchy boundary condi- 
tions, one has nothing further to do (unless one is studying an open system, in which 
case there remains the integration over the initial data of the environment). And since 
the classical trajectories of a system are easy to associate with a consistent physical 
intuition, one can not only compute the correct final answer, but also appreciate the 
intermediate processes that generate the ultimate result. Since most linear problems 
are essentially of conceptual interest, this is certainly a consummation to be wished: 
Classical mechanics can greatly simplify quantum mechanics. 

In the linear regime, the converse can also be true. It is widely assumed that, 
even with the simple environmental model of independent oscillators, the classical 
Fokker-Planck equation for a Brownian oscillator is only local in time if one assumes 
a particular spectral density for the environmental bath. This is not true. By ap- 
plying an insight from a recent solution of the quantum problem, we can abandon 
the Markovian assumption in classical Brownian motion, and still obtain a local 
Fokker-Planck equation. In the non-Markovian regime, the correct identification of 
the physical noise and the effective application of the fluctuation-dissipation relation 
become much more complicated || . Without the use of methods inspired by quantum 
mechanics, it would be very difficult indeed to discern the Fokker-Planck equation: 
In this problem, quantum mechanics can greatly simplify classical mechanics. More- 
over, the demonstration that the classical problem is indeed dynamically identical to 
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the quantum problem must remove any lingering suspicion that influence functionals 
(or any other essentially quantum techniques) are necessary to capture the quantum 
dynamics. 

To proceed to explicit details, we recall the Hamiltonian for the Brownian har- 
monic oscillator linearly coupled to a bath of independent oscillators: 
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where g u is the coupling constant, whose u;-dependence effectively incorporates the 
spectral density of the environment. We also recall the Liouville propagator for the 
probability density in phase space: 

f(Q F ,P F ; {q^ F , PujF };t) = JdQjdPjDqjDpjSiPp - P{t))5{Q F - Q{t)) 

X / 11%^ ~ <lw(f))f>(PuF -Pw(t)) 

xf{Qi,Pi] {q u i,Pui},ti) , (2) 

where Q(t), qu(t) and P(t),p tJ (t) denote the values of the canonical variables at time 
t after evolution, under Hamilton's equations, from the initial values Qi, Pi and 
{qui,Pwi} at time £/. We write Q(t) instead of Q(t; Qi, Pi, {qui,Pui}), etc., only in 
order to avoid an unwieldy notation. It will be important to remember, despite this 
shorthand convention, that Q(t) and P(t) are time-dependent linear combinations of 
all the initial variables. 

As in the standard theoretical version of Brownian motion, we assume that the 
intitial ensemble factorizes and that its environmental part is Gaussian^ 

f{Qi,Pi; {qi{u), P i{u)},ti) = f(Qi,Pi;ti) n(^)e-^ 2 ^ + ^l . ( 3 ) 

Since we wish to know only the reduced distribution f(Q,P;t) over the observed 
sector of phase space, we integrate over the environmental sector in (Q). With the 
factorized initial condition (|3|), this is equivalent to propagating / using the equations 
of motion P(t) = mQ{t) and 

Q(t) + n 2 Q(t) + K(t)Qi + / dt'K(t - t')Q(t') = -A2 , (4) 

Jo m 

where K(t) = (1/m) J Q °°duj g 2 cos cut, and F(t) is a stochastic force. Since the environ- 
mental part of @ is Gaussian, F is completely described by its two-point correlation 

1 Slight variations in this initial distribution may be contemplated, as may variations in the 
Hamiltonian that are obtained by canonical transformations mixing the system and environmental 
sectors. These variations can be encompassed by straightforward generalizations of the present 
discussion. 



function: 

(F(t)F(t')) = [°°du ^ coscu(t - t') . (5) 
Jo p u 

Eqn. (Q) is the (generalized) Langevin equation. The standard classical approach 
to this problem relies on separating the systematic and stochastic terms in this equa- 
tion 0]. In the special case where ng^ = 4m7 is constant, the integral kernel K be- 
comes a delta function, and one obtains a stochastic differential equation with Ohmic 
damping. The entire LHS of (f|) is then systematic, and the RHS can be unambigu- 
ously identified as the noise. The usual classical prescription (3 U = (fcgT) -1 , indepen- 
dent of to as well, also ensures that the environmental noise is white. In this case, 
it is easy to use the fluctuation-dissipation relation (F(t)F(t')) = Am^kBT5{t — t') 
to derive an evolution equation for the reduced phase space distribution /, which is 
local in time: the Fokker-Planck equation. 

For general g w , however, Langevin equations are integro-differential, incorporating 
a back-reaction term with memory. The noise is also colored in general {i.e., the 
correlation function is not a delta function), even with constant (3 U . One might 
think that these circumstances would prevent the generalized Fokker-Planck equation 
from being local. In fact, because the memory term depends on the past history of 
Q, over all of which it has been affected by the environmental noise, it actually 
contains an implicit stochastic part. So the physical noise receives a contribution 
from back- reaction, as well as the explicit "bare noise." The crucial identification of 
the systematic and stochastic parts of the Langevin equation is therefore more subtle. 

One can still derive a generalized Fokker-Planck equation by the standard method, 
using the colored physical noise that includes back-reaction, and the non-local sys- 
tematic evolution. Remarkably, it will turn out that a rather complicated kind of 
fluctuation-dissipation relation can be found, which will guarantee that the general- 
ized Fokker-Planck equation derived from Eqn. (J|) will always be local in time ||. 
But this approach is quite cumbersome. 

A much shorter derivation is available, inspired by the quantum mechanical treat- 
ment of Ref. 0]. Integrating over the environmental sector in Eqn. (^), and using a 
representation for the two surviving delta functions that is more familiar in quantum 
contexts, we can write 

f{Q F ,P F ;t)=i JdQjdPjfiQj^^tj) jDd qi D Pl e -U~<^«x+^] 

(6) 



Because the Langevin equation is linear, the solutions Q(t),P(t) depend linearly on 
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the initial variables {qw,£w} and Qj,Pj. 

The k and k! pre-factors in @ can be replaced by derivatives under the integral 
with respect to Qp and Pp. The thus- restored delta functions can then be used to 
replace the Qi and Pj terms in Q(t) and P(t) with linear combinations of Qp, Pp, and 
the and p^i- The Gaussian integrals over the environmental variables can then 
be performed. The prefactors proportional to q u j and p^i lead to terms proportional 
to k and k', which can once again be expressed as derivatives with respect to Qp and 
Pp. At the end of this process, we drop the subscripts on Qp and Pp, and obtain an 
equation of the form 

The time-dependent co-efficients may be expressed in terms of the parameters in the 
Hamiltonian (p. Their calculation is straightforward but tedious, and we will not 
perform it here. Instead we merely note that time-dependent co-efficients, as well 
as the so-called "anomalous diffusion" term proportional to d(t), both previously 
encountered in quantum treatments, do also appear in the classical problem. 

If we set P n = ^-tanh^^ in @ and @, then equation (^) is the Hu-Paz- 
Zhang master equation J7J, in the Wigner representation. Any positive choice of 
/3 n , however, will still provide a local equation for/. And while (3 n = (ksT)" 1 is 
prescribed by classical statistical mechanics (and for Ohmic spectra this differs from 
quantum mechanics in permitting white noise for all temperatures), any positive (3 n 
still defines an allowed (i.e., positive) classical ensemble density for the environment. 
In the case of factorized, environmentally Gaussian initial conditions, then, classical 
and quantum Brownian motion can be said to differ only in the constraints prescribed 
for f(Qi,Pi;ti). Dynamically, the quantum and classical cases are identical. 

Moving beyond the specific linear problem of Brownian motion, it is straightfor- 
ward to verify that, for the most general time-dependent quadratic Hamiltonian, the 
propagator for the quantum mechanical Wigner function is of exactly the same form 
as Eqn. d). This means that the quantum Liouville equation for a linear system in 
the Wigner representation will be the same as the classical Liouville equation for the 
corresponding classical Hamiltonian. Moreover, this means that the classical equa- 
tions of motion, with their classically causal boundary conditions, provide as exact 
and complete an expression of linear quantum dynamics as any operator equation 
or path integral. And for linear open quantum systems with an initially Gaussian 
environment, the reduced dynamics is directly expressed in the c- number Langevin 
equations. 

These c-number equations are both conceptually and computationally much less 
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demanding than the more general formulations of quantum dynamics. If one wishes to 
obtain the final Wigner function f(Q, P; t) of a linear quantum system, one need only 
solve the time-reversed classical equations of motion to determine the initial phase- 
space point corresponding to the final point (Q,P). The initial Wigner function at 
this initial point is the quantity sought. The Green's functions needed to do this for 
final Cauchy conditions are typically much easier to use than those required for the 
mixed initial and final conditions of other approaches, and since the initial Wigner 
function is given, one need do no further work once the classical equations are solved. 
And of course this program remains effective for linear systems with any number of 
degrees of freedom. 

Furthermore, in many cases one may be able to gain considerable information from 
the equations of motion alone, with only qualitative discussion of the initial Wigner 
function. For example, knowing that a bath of independent oscillators provides a 
Langevin noise, one can conclude that environmental fluctuations will tend to even 
out oscillations in the Wigner function of a Brownian oscillator. Rapid oscillations 
between positive and negative values will be suppressed quickly, and this provides one 
of the simplest intuitive pictures of the process of decoherence (since such oscillations 
are the Wigner signature of interference between classically distinct pointer states). 

As a second and less trivial example, we can also use classical Langevin equa- 
tions to understand why so-called "supra-Ohmic" environments tend to induce little 
decoherence in Brownian oscillators || (a more complete discussion will be given else- 
where II). In the present formalism, supra-Ohmic environments are those for which 

increases with uo (up to some cut-off scale, generally taken to be much larger than 
all other scales in the problem). If one compares Ohmic and supra-Ohmic models 
in which are equal, one typically finds that (after a few cut-off timescales) the 
supra-Ohmic environment produces weaker diffusion in the Brownian oscillator. This 
is because stronger coupling at high frequencies means that adiabatic dragging of fast 
environmental oscillators gives the Brownian particle a higher effective mass. Indeed, 
Brownian motion in supra-Ohmic environments heavily weighted in the ultraviolet 
rapidly becomes indistinguishable from Ohmic motion with a renormalized mass || . 

One can easily anticipate this effect by using an adiabatic approximation to write 
an effective Lagrangian without the fast degrees of freedom, but it is interesting to see 
how it appears in the exact solution of the full problem. Computing the co-efficients 
in the master equation (^) merely confirms that diffusion dies away; studying the 
Langevin equation @ reveals how this occurs. After solving the equation, one can 
always in principle re-write it in the form 

Q(t) + n 2 (t)Q(t) + 2j(t)Q(t) = i[F(t) + F BR (t)} , (8) 
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where in addition to the bare force F (whether stochastic or systematic) there also 
appears the back-reaction force Fbr- In the case of a UV-dominated supra-Ohmic 
environment, there is a brief inertial epoch that lasts for a few cut-off times, during 
which the system behaves much as an Ohmic model, and the back-reaction force is 
small. After this epoch, however, a remarkable regime emerges in which any impulse 
applied to the Brownian oscillator as a bare force is rapidly echoed in the back-reaction 
force, with opposite sign and nearly equal magnitude ||. This "counter-punch" ef- 
fectively suppresses all forces acting on the oscillator, unless they vary rapidly on the 
cut-off timescale, by a large factor equal to the ratio of bare and renormalized masses 
found in the adiabatic analysis. Thus we achieve a nice intuitive picture, in which 
the supra-Ohmic environment appears to deliberately isolate the Brownian particle, 
and in which the importance of back-reaction is clearly exhibited. 

When using classical intuition in quantum problems, one must of course be careful. 
Wigner functions cannot, in general, be interpreted as probability densities, because 
they can be negative. Even when one happens to have a positive definite Wigner 
function, as initially one has in the environmental sector for most Brownian motion 
problems, one cannot assume that the classical interpretation will persist for all times. 
For example, the picture just presented of decoherence, as due to classically stochas- 
tic environmental noise, may seem to disagree with the usual explanation in purely 
quantum language, according to which the different branches of a "Schrodinger's Cat 
state" excite orthogonal states of the environment ||. 

Actually, the two explanations can be reconciled by realizing that our classical 
description effectively used a mixture of the Schrodinger picture for the Brownian 
system and the Heisenberg picture for the environment. It was thus able to refer 
only to the initial Wigner function of the environment, which was interpretable as a 
probability density. But if we were to use a uniformly Schrodinger picture, and follow 
the evolution of the entire Wigner function before integrating over the environmental 
sector, we would see that the negative values in the initial Wigner function of the sys- 
tem spread rapidly into the environmental sector. When the environmental integrals 
are performed, these contaminating oscillations into negative values provide cancel- 
lations that are the exact analogue of the vanishing inner products in the standard 
discussion. So the fact that the environment can affect the observed system only as a 
stochastic force does not allow us to conclude that the environment actually remains 
in a state described by a probabilistic ensemble. 

Even if a Wigner function is positive definite, it still obeys a non-classical con- 
straint associated with the Uncertainty Principle: 




(9) 
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where / dT denotes the integral over all of phase space. This constraint has important 
consequences. For example, a positive definite Wigner function has a well-defined 
Boltzman entropy, but this entropy will generally be larger than the von Neumann 
entropy. This can be interpreted as being due to the fact that nearby points in phase 
space are not distinct events in quantum mechanics, but are in a sense parts of the 
same object. So the von Neumann entropy represents a special kind of coarse-grained 
entropy. 

(This idea provides a way to express the distinct quantum mechanical concepts of 
correlation and entanglement, in the Wigner representation. Compare the following 
two Wigner functions, each for a system of two oscillators: 

+ e -^(Ql+a) 2 e -^(Qi+a) 2 ] 

f 2 = z 2 e-^ p ? +p 2 ) e- n ir( Q i +Q 2 ) e rQlQ2 . (10) 

Classically, one would say that both of these functions describe correlations between 
oscillators 1 and 2. But quantum mechanically, the first describes a correlation (for 
a 2 >> while the second describes an entanglement. We can conclude that 

entanglement is a correlation between parts of a single quantum mechanical object 
T0|| ! an d that the distinction between correlation and entanglement is not a separate 



piece of "quantum magic" from the uncertainty relation, but is actually a consequence 
of it.) 

Despite the subtleties involved in interpreting Wigner functions, we emphasize 
that as long as one only uses classical intuition to understand how the classical equa- 
tions of motion propagate them from initial to final times, there is nothing that can 
go wrong. This careful use of classical intuition can be of considerable value in appre- 
ciating the physics of linear quantum systems. One can expect to shed new light on 
the meaning of some important quantum mechanical effects, by identifying aspects of 
classical dynamics that are responsible for their generation. This is the good news, 
but it has a more pessimistic complement. 

Because illustrations in linear models of basic quantum concepts are usually 
treated using techniques applicable to non-linear models as well, it is easy to re- 
ceive an impression that these models are providing an insight into essentially quan- 
tum physics. It is important to realize, however, that the only specifically quantum 
aspects of any of these models reside in the initial conditions. To probe genuine quan- 
tum dynamics, we must accept the challenge of non-linearity, and must not confuse 
generality of formalism with generality of physics. 
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